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Abstract 

Cross sections and photon asymmetries of the exclusive 16 0(7, pn) 14 N and 16 0(7, pp) 14 C knock- 
out reactions are calculated for transitions to the low-lying discrete final states of the residual 
nucleus in the photon-energy range between 100 and 400 MeV. Exclusive reactions may represent 
a test of reaction mechanisms and a promising tool for investigating the dynamics of nucleon pairs 
^ \ in different states. Cross sections and asymmetries for both (7,pn) and (7,pp) turn out to be only 

slightly affected by short-range correlations and dominated by two-body currents. Therefore, two- 
nucleon knockout reactions induced by real photons appear well suited to investigate the nuclear 
current and the selectivity of individual transitions to its different components. 
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g . 1. INTRODUCTION 

> : 

^ ■ Electromagnetically induced two-nucleon knockout reactions have long been considered 

an important tool for investigating the properties of nucleon pairs within nuclei and their 
interaction at short distance |j],[2|]. 

Direct information on dynamical short-range correlations (SRC), which are linked to the 
short-ranged repulsive core of the NN-interaction, can be obtained if one assumes that the 
real or virtual photon hits, through a one-body current, either nucleon of a correlated pair 
and both nucleons are then ejected. However, nucleon pairs can also be ejected by two- 
nucleon currents, which effectively take into account the influence of subnuclear degrees of 
freedom like mesons and isobars. A reliable treatment of both these competing processes 
as well as of other reaction mechanisms is in principle needed before one can draw definite 
conclusions in comparison with data. However, since the role and relevance of these contribu- 
tions is different in different reactions and kinematics, it is possible to envisage appropriate 
situations where various specific effects can be disentangled and separately investigated. 

The most promising tool for studying SRC is represented by the (e,e'pp) reaction, where 
the effect of two-body currents is less dominant as compared to the (e,e'pn) and (7,NN) 
processes. The first exploratory triple coincidence measurements of the (e,e'pp) cross section, 
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performed at NIKHEF on 12 C |3],|4j], gave clear evidence of two-nucleon knockout in the dip 
region. Peaks corresponding to knockout of protons from the lp and Is shells were clearly 
recognized in the missing-energy spectrum. The measured cross sections were found in 
satisfactory agreement with the model calculations of refs. [||,|| in the region corresponding 
to (lp) 2 knockout and indicated that a large part of the strength can be attributed to SRC. 
In most recent experiments on 16 O at NIKHEF |7||| and Mainz || it was possible to achieve 
sufficient energy resolution to allow the separation of the cross section related to distinct 
states of 14 C. A further experiment on 16 with improved statistics has been approved 
10]. The presence of discrete final states with well-defined angular momentum 

since specific final states may act as filter for 



in Mainz 

makes 16 O a particularly attractive target 
the study of different reaction processes involving SRC and two-nucleon currents |]TT,12 
Moreover, the first and so far only available calculation of the two-nucleon spectral function 

16Q 



for a finite nucleus, has been applied just to 
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The corresponding pair removal 
amplitudes have been included in the description of the (e,e'pp) knockout reaction ||T2"f . 
The predicted selectivity of the considered reaction involving different final states, which is 
confirmed by the fair agreement obtained in the first comparisons with data |||9|], opens up 
good perspectives for the study of SRC in (e,e'pp) reactions. 

The (7,NN) reactions appear dominated by medium-range meson-exchange currents and 
a good understanding of the dominant single-pion-exchange mechanisms is essential before 
one can hope of investigating shorter range effects [0,0 • 

First exploratory (7,NN) experiments performed in Bonn |16j and Tokio [17| had neither 
good statistical accuracy nor sufficient energy resolution to allow a straightforward inves- 
tigation of direct processes. Early work in Mainz ]18[ had good energy resolution allowing 
the identification of nucleons emitted from different shells. Further experiments with good 
statistics and resolution were carried out in Mainz |T^|,^DJ and Lund |?T| in order to disen- 
tangle the reaction mechanisms. A major conclusion was obtained through comparison of 
missing-energy spectra with the results of the theoretical model of ref. [2~2~|, where different 
contributing mechanisms are included in the description of the pion and nucleon photopro- 
duction channels in a unified diagrammatic approach: only the low missing-energy region 
(E 2m < 40 MeV) is predominantly fed by direct photoabsorption on a nucleon pair, whereas 
for E 2m > 70 MeV other mechanisms, e.g. contributions originating from quasifree (7,71") 
production and subsequent pion final-state interactions, become dominant. This result im- 
plies that only low values of E 2m must be considered when studying the properties of nucleon 
pairs and was qualitatively confirmed by the analysis of the data in ref. [2~0| through the 
model of ref. |14| . 

Measurements with sufficient energy resolution to separate the low-lying final states 
of the residual nucleus are of great interest also in the case of photon-induced reactions. 
They represent a stringent test of reaction models and allow one to investigate the state 
dependence of the photoabsorption process on nucleon pairs. However, most of the data 
collected up till now could not separate individual final states and simultaneously provide 
dependences of the cross sections on dynamical variables. The high-resolution experiments 
performed in Lund [^T| for the 16 0(7, pn) 14 N reaction could separate final states, but, owing 
to low statistics, no momentum distribution for even the strongest states could be obtained. 
A new high resolution experiment for the 16 0(7,pn) 14 N and 16 0(7,pp) 14 C reactions in the 
photon-energy range between 90 and 270 MeV, aiming at separating final states and at 
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determining also the momentum distributions for each state, has been recently approved in 
Mainz |?3f| . 

Various theoretical models have been developed over the last years to deal with two- 
nucleon emission processes. The model of ref. |22| aims at the best microscopic description 
of the photoabsorption process and includes many different reaction mechanisms. This ap- 
proach is extremely useful to understand and disentangle the main mechanisms contributing 
in different kinematical regions and is able to reproduce their general features. However, 
it cannot describe the details of the nuclear structure aspects of the reaction, which are 
dealt with in a nuclear-matter approach with a local-density approximation. A detailed 
description of these aspects is essential if one wants to investigate the conditions of pairs of 
nucleons in finite nuclei and their state dependence through exclusive knockout reactions. A 
more proper description of nuclear structure properties is given by the two models developed 
in refs. ]5],|| and |14]]T5|| , where, on the other hand, only the direct two-nucleon knockout 
mechanism is considered. Both of them include SRC, one-body and two-body currents and 
final-state interactions in a completely unfactorized calculation. These different aspects have 
been investigated at various steps and with somewhat different theoretical ingredients in the 
two models, which have both been applied to the analysis of existing data for electron- and 
photon-induced reactions. 

In refs. jnj.15] reactions to discrete final states are not explicitly investigated but in 
very particular situations. In refs. |||§ the various approximations used in the calculations 
make it difficult to evaluate the cross sections for transitions to individual final states. 
Thus a sum over all contributions given by different pairs of nucleons in the same shell is 
performed, which experimentally corresponds to detect all the pairs of nucleons coming out 
of the considered shell. This assumption is largely compatible with the energy resolution 
of the first experiments, but is not adequate for the analysis of the most recent exclusive 
data already available for the (e,e'pp) reaction |7|-|| and of the new data that are becoming 
available for both electron- and photon- induced reactions pH|J2^ |. 



Exclusive (e,e'pp) reactions have been investigated in refs. |pT,l2[. In ref. JTT| the the- 
oretical model of refs. || has been improved with a full treatment of antisymmetrization 
and of spin and isospin couplings. In ref. [|T^] the pair removal amplitudes obtained from 
the calculation of the two-proton spectral function of 16 O have been included in the 
reaction calculation. The treatment of the A isobar current is more specifically discussed in 
ref. ||24j| , where some first numerical results for the exclusive (7,pp) reaction are given but 
the spectral function is not included in the model. 

In this paper results for the exclusive 16 0(7,pn) 14 N and 16 0(7,pp) 14 C reactions to in- 
dividual low-lying discrete final states of the residual nucleus are presented and discussed. 
These two reactions appear of particular interest. Exclusive data for both of them are 
expected from an approved experiment in Mainz |23| . 

The approach already applied to the 16 0(e, e'pp) 14 C reaction |12[ is here applied to the 
16 0(7, pp) 14 C reaction. In this approach the pair removal amplitudes are obtained from the 
calculation of the spectral function, where both long-range and short-range correlations are 
treated consistently. Different relative and center-of-mass (CM) states of the pair contribute 
to a specific transition. They can be separated in the calculation and their effects can 
individually be investigated. A calculation of the spectral function for a pn pair is not 
available yet. Thus for the 16 0(7,pn) 14 N reaction the simpler prescription of ref. |TT[] is 
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adopted here, where the two-nucleon overlap is given by the product of a coupled and 
fully antisymmetrized pair function of the shell model (SM) and a Jastrow type correlation 
function. In the two-body current for a pn pair contributions given by both seagull and 
pion-in-flight diagrams have been added to the A isobar current. This represents a further 
improvement with respect to the previous calculations for (7,pn) reactions of ref . |6j , where 
pion-in-flight diagrams were neglected. With respect to ref. ||15|| , which is mainly devoted 
to investigate the quasi-deuteron kinematics and the behaviour of nucleon pairs knocked- 
out from different shells, this paper aims at investigating the transitions to discrete states, 
in order to explore their selectivity to the different components of the nuclear current. 
Moreover, the comparison of the results obtained with different models can give a deeper 
insight into the mechanisms of two-nucleon emission and the nuclear structure properties. 

In this paper numerical results of cross sections and asymmetries are presented in the 
photon-energy range between 100 and 400 MeV. It was already observed in ref. || how 
the asymmetry of the cross section, that can be measured with linearly polarized photons, 
is helpful to investigate the reaction mechanism and the role of the different terms of the 
nuclear current. First asymmetry measurements for the 16 0(7,pp) and 16 0(7,pn) reactions 
have been carried out at LEGS |25] in the photon energy range between 245 and 315 MeV. 
At these energies, around the peak of A resonance, the contribution of the isobar current is 
dominant and the calculations are extremely sensitive to its treatment. The data indicate the 
need of a careful description of this important theoretical ingredient. In this first experiment, 
however, the energy resolution was not sufficient to separate final nuclear states. 

The more refined treatment of spin degrees of freedom in the present approach is expected 
to give important effects on a polarization observable such as the asymmetry. It seems 
moreover interesting to investigate the dependence of this observable on the different states 
of the nucleon pair. 

The theoretical approach is outlined in sect. 2. Numerical results for the 16 0(7,pn) 14 N 
and 16 0(7,pp) 14 C reactions are presented and discussed in sect. 3 and sect. 4, respectively. 
Some conclusions are drawn in sect. 5. 



2. THEORETICAL APPROACH 



The coincidence cross section for the reaction induced by a photon, with energy E 1 , 
where two nucleons, with momenta p[, and p 2 and energies E[ and E f 2 , are ejected from a 
nucleus, is given, after integrating over E' 2 , by 

^ = ^O f /-^ T (1) 

dEidQW 2E„ ec {) 

where f2f = •p' x E' x j)'< 1 E'< 1 is the phase-space factor and integration over E' 2 produces the recoil 
factor 



f- 1 — 1 _ El. p ' 2 ' Pb (o\ 
/rec ~ E B |p' 2 |2 ' 
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where Eb and p B are the energy and momentum of the residual nucleus. The transverse 
structure function W^, which only depends on E 7 , p[, p' 2 , the angles 71, between the mo- 
mentum of the incident photon q and p[, 72, between q and p' 2 , and 712, between p[ and 
p' 2 i is expressed in terms of the components of the hadron tensor W^ u ||, i.e. 

W T = W xx + W yy (3) 

and is thus given by bilinear combinations of the Fourier transforms of the transition matrix 
elements of the nuclear current density operator taken between initial and final nuclear 
states. 

If the residual nucleus is left in a discrete eigenstate of its Hamiltonian, i.e. for an 
exclusive process, and under the assumption of a direct knockout mechanism, the transition 
matrix elements can be written as |o.]J 



J(q) = J ^i{r 1 a- 1 ,r 2 (r 2 )J(r,r 1 (T 1 ,r2(T2)^i{r 1 (T 1 ,r2(T2) 

x e l9 ' r drdridr 2 dcr 1 dcr2. (4) 

Eq. dj) contains three main ingredients: the two-nucleon overlap integral the nuclear 
current J and the final-state wave function ipf. 

The derivation of Eq. (f|) involves bound and scattering states, ipi and if)f, which are 
consistently derived from an energy- dependent non-hermitian Feshbach-type Hamiltonian 
for the considered final state of the residual nucleus. They are eigenf unctions of this Hamil- 
tonian at negative and positive energy eigenvalues, respectively In practice, it is not 
possible to achieve this consistency and the treatment of initial and final state correlations 
proceeds separately with different approximations. 

In the final-state wave function ipf each of the outgoing nucleons interacts with the resid- 
ual nucleus while the mutual interaction between the two outgoing nucleons is neglected. 
The scattering state is thus written as the product of two uncoupled single-particle distorted 
wave functions, eigenfunctions of a complex phenomenological optical potential which con- 
tains a central, a Coulomb and a spin-orbit term. 

The two-nucleon overlap integral ipi contains the information on nuclear structure and 
allows one to write the cross section in terms of the two- hole spectral function [Q]. Since 
only a calculation of the two-proton spectral function of 16 is available at present [13 , 12| , 
two different prescriptions have been here adopted for pp and pn knockout. 

For the 16 0(7, pp) 14 C reaction the two-nucleon overlap integrals are taken from the calcu- 
lation of the spectral function [jT3,|T3] • They have been obtained from a two-step procedure, 



where long-range and short-range correlations are treated in a separate but consistent way. 
The calculation of long-range correlations is performed in a SM space large enough to in- 
corporate the corresponding collective features which influence the pair removal amplitudes. 
The single-particle propagators used for this dressed Random Phase Approximation descrip- 
tion of the two-particle propagator also include the effect of both long-range and short-range 
correlations. In the second step that part of the pair removal amplitudes which describes 
the relative motion of the pair is supplemented by defect functions obtained from the same 
G-matrix which is also used as the effective interaction in the RPA calculation. 

For a discrete final state of the 14 C nucleus, with angular momentum quantum numbers 
JM, the two-nucleon overlap integral is expressed in terms of relative and CM wave functions 
as 
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^(ri<ri, r 2 <T 2 ) = J2 c idSjNL<t>nlSj(ri2)RNL(R) [S/s( fi r,CTi,(T 2 )y L (fi R )] JM , (5) 

nlSjNL 

where 

D ( r l+ r 2) ,„x 

r 12 = r 1 -r 2 , = (6) 

are the relative and CM variables. The brackets in Eq. @ indicate angular momentum 
coupling of the angular and spin wave function 3 of relative motion with the spherical 
harmonic of the CM coordinate to the total angular momentum quantum numbers JM. 
The CM radial wave function Rnl is that of a harmonic oscillator (h.o.), with oscillator 
parameter b = 1.77 fm. SRC are included in the radial wave function cf) of relative motion 
through a defect function defined by the difference between <ft and the uncorrelated relative 



h.o. wave function of the pair R nh i.e. [ fL3[] 



<t>niSj{ri2) = Rni(r 12 ) + D lSj {r l2 ). (7) 
These defect wave functions were obtained by solving the Bethe-Goldstone equation in mo- 



mentum space for 16 [26]. They depend on different quantum numbers I, S, I. In the 
following the partial wave notation 2S+1 lj, for I = S, P, D, is used for the relative states. 

The coefficients c 1 in Eq. fl5|) contain contributions from a SM space which includes the 
Os up to the lpOf shells. More details are given in ref. [13 and in ref. [I2|, where the same 



approach is applied to the reaction 16 0(e, e'pp) 14 C. 

For the reaction 16 0(7, pn) 14 N, where a calculation of the spectral function is not avail- 



able, the simpler prescription already applied in ref. [11| to (e,e'pp) has been adopted here. 
The two-nucleon overlap is thus given by 

ipi(ri(T U r 2 (T 2 ) ~ $jA/(r'icri,r2cr 2 )/(r 1 2)X TT 3(r 1 ,r2), (8) 

where $jm^tt 3 is the coupled and fully antisymmetrized pair function of the shell model, 
X TT3 is its isospin part and / is a correlation function of Jastrow type which incorporates 
SRC. Only the central part of the correlation function is retained in the calculations. Thus 
long-range correlations are not included in this simpler approach. The overlap function of 
Eq. (^j) describes a pair of nucleons ejected from the same or different shells. The model 
includes the spectral strength A. In the calculations we have assumed A = 1. 

The nuclear current operator in Eq. (|]) is the sum of a one-body and a two-body part. In 
the one-body part convective and spin currents are included. The two-body current is derived 
from the effective Lagrangian of ref. [^7[, performing a non relativistic reduction of the lowest- 
order Feynman diagrams with one-pion exchange. We have thus currents corresponding 
to the seagull and pion-in-flight diagrams and to the diagrams with intermediate isobar 
configurations, i.e. 

J (2) (r, n<ri, r 2 <r 2 ) = J sca (r, r x a x , r 2 a 2 ) + J w (r, r x a x , r 2 a 2 ) 

+ J A (r,ri<ri,r a o-3)- (9) 

In the coordinate space the seagull and pion-in-flight currents are 
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J^(r,r l( r u r 2( T 2 ) = - f — (n x r 2 ) 3 [o-«<f(ri - r){a^ ■ r 12 )\ 

/ 1 \ e _Atri2 

x 1 + + (l<-2), (10) 

V F12/ fJ>rn 
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J*(r, r 1 cr 1 , r 2 er 2 ) = ~J^2 ( Tl x Vl ( <Tl ' Vl ) (°" 2 ' V2 ) 

p -/x|ri-r| p -M|r 2 -r| 

r+ (1-2), (11) 

yu|rx — r\ fj,\r 2 — r\ 

where / 2 /(47r) = 0.079 and \i is the pion mass. 

The operator form of the A current has been derived in ref. [21|. It is given by the 



sum of the contributions of two types of processes, corresponding to the A-excitation and 
A-deexcitation currents. The first process (I) describes A-excitation by photon absorption 
and subsequent deexcitation by pion exchange, while the second (II) describes the time 
interchange of the two steps, i.e., first excitation of a virtual A by pion exchange in a 
NN collision and subsequent deexcitation by photon absorption. The propagator of the 
resonance, G A , depends on the invariant energy a/s of the A, which is different for processes 
I and II. For the deexcitation current the static approximation can be applied, i.e. 

G£=(M A -M)-\ (12) 



where M A = 1232 MeV. For the excitation current we use 28 



Gi = (m a - v^I - ^r A (v*)) 1 , (13) 

with 

/ si= y /s^-M, (14) 



where ^/snn is the experimentally measured invariant energy of the two outgoing nucleons 
and the energy- dependent decay width of the A, r A , has been taken in the calculations 
according to the parameterization of ref. . 
The sum of the two processes gives 

J A (r, r x <r x , r 2 a- 2 ) =j6(r- n) {i (G& + G A ) [4r 2i3 A(r 12 , cr 1; <x 2 ) - (r x x r 2 ) 3 
x B(r 12 , (T U (T 2 )\ + 2{G\ - Gl) [(n x r 2 ) 3 A(r 12 , a x , a 2 ) 

+ T2,3%2,(T 1 ,(r 2 )]} + (1-2), (15) 

where 

A(r 12 , <t i , <x 2 ) = (q x f 12 ) (<r 2 ■ f 12 ) r {1) (n 2 ) - (g x <x 2 ) F (2) (r 12 ), (16) 

B(ri2, <ri, a- 2 ) = qx (o-i x f 12 ) (<r 2 • f 12 ) F (1) (ri 2 ) - g x (tn x cr 2 ) F (2) (ri 2 ), (17) 
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\ //ri 2 l~L z r 12 z J nr 12 

— + ^—2 " > ( 19 ) 

and the factor 7 collects various coupling constants 

_ f^/NA fnNN fixNA 
7 ~ 367T/i ' 1 J 

The charge-exchange term (ri x T2)3 in the two-body current vanishes for a pp pair. 
Thus only a part of the A current in Eq. contributes to (7,pp) reactions, whereas all 
the terms in eqs. (jT0,|rTL[T5[) contribute to (7,pn). On this basis it is generally expected that 
the two-body current is much more important for the knockout of a pn pair than of a pp 
pair. 

From the analysis of the isospin matrix structure of processes I and II, it can be under- 
stood JHJ] that for photon absorption on a pn pair only the excitation current contributes in 
a state with T = 0, whereas only the deexcitation current contributes in a state with T = 1. 
Since the propagator gives rise to a pronounced resonant behaviour of the matrix ele- 
ments of the excitation current in the photon-energy region between 200 and 400 MeV, we 
expect that in this region the contribution of the A current for the knockout of a pn pair 
with T = is much larger than for a pn pair with T — 1. 

Both processes I and II contribute to photon absorption on a pp pair (T = 1). For 
the knockout of a pp pair in a relative 1 Sq state the generally dominant magnetic dipole 
NN <-> iVA transition is suppressed because of total angular momentum and parity conser- 



vation |30j . In our calculations contributions of all multipoles are included through the term 
c i( l r , which represents the incident photon in the transition amplitudes of Eq. (H). Only if 
we set e i( l r = 1 the calculation for the A current is restricted to a magnetic-dipole transition. 
In any case the effect of higher multipoles is generally smaller and for a 1 Sq pp pair the A 
current, although nonvanishing, turns out to be much less important than for other relative 
states. This was already observed in ref. [|1^] for the (e,e'pp) reaction. In that case the 
contribution of the excitation current, which is generally dominant, is strongly reduced on a 
1 S'o pp pair, where it becomes of about the same size or even smaller than that of the deex- 
citation current, which is generally small. As a consequence, in (e,e'pp) reactions situations 
where the removal of a 1 Sq pair is dominant are particularly well suited to emphasize and 
thus investigate effects of SRC |7|,[12|. This conclusion cannot be directly applied to (7,pp) 
reactions, where only the transverse component of the nuclear response contributes and the 
longitudinal component, where the effects od SRC show up more strongly, is absent. 

In photon-induced reactions with linearly polarized photons one can measure the photon 
asymmetry S, given by the difference between the cross sections with linear photon polar- 
ization parallel and perpendicular to the reaction plane. It can also be expressed in terms 
of the ratio between the two transverse structure functions Wtt and Wt Pfl, as 

= w xx - wyy = Wtt , . 

w xx + wyy w T ' 

The asymmetry is particularly sensitive to spin variables and thus to the different terms 
of the nuclear current and to their interference, whose effects on Wt and Wtt are emphasized 
in the ratio, where, on the contrary, other effects are smoothed out |6|jr5|]. 
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In the following numerical results are presented of cross sections and asymmetries for 
the 16 0(7,pn) 14 N and 16 0(7,pp) 14 C reactions for transitions to discrete final states of the 
residual nucleus. High-resolution experiments are needed to separate individual final states 
in the excitation-energy spectrum. In our calculations each state is characterized by a 
particular value of the missing energy, given by 

E2m = E 1 — T[ — T' 2 — T B = E s + E x , (22) 

where T{, T' 2 and Tb are the kinetic energies of the two outgoing nucleons and of the residual 
nucleus, respectively, E s is the separation energy at threshold for two-nucleon emission and 
E x is the excitation energy of the residual nucleus. 



3. THE REACTION 16 0(7,pn) 14 N 



The theoretical approach for two-nucleon knockout outlined in sect. 2 has been applied 
to the reaction 16 0(7, pn) 14 N, for transitions to discrete final states of the residual nucleus. 
Only the states with an excitation energy lower than 12 MeV have been considered in the 
calculations, according to the experimental analysis of ref . |]19[ , where the direct two-nucleon 
knockout mechanism is well established for missing energies up to 40 MeV in 12 C. Therefore 
we consider here transitions to the if (T = 0) ground state of 14 N, to the + (T = 1) 
state at 2.31 MeV, to the (T = 0) state at 3.95 MeV, to the 2+ (T = 0) state at 7.03 
MeV and to the 3 + (T = 0) state at 11.05 MeV. We do not consider the negative parity 
states at low energy and other positive parity states, as they are not assumed as two-hole 
states in the independent particle model and therefore are not likely to be excited in a direct 
knockout [ |31|j . 

The two-nucleon overlap is taken in the calculations according to the prescription of 
Eq. (||). In the SM wave function, all the considered states are described as two- hole states 
in the p shell: 1+ and + are (pi/2)~ 2 holes, 1 2 and 2 + are (pi/2P3/2) _1 holes and 3 + is a 
{p3/2)~ 2 hole. In the frame of a h.o. basis, their different components in the relative and 
CM motion are given in Table I with the corresponding amplitudes. This decomposition 
can be helpful to understand the general features of cross sections and asymmetries. In 
the calculations, however, the more realistic Woods-Saxon single-particle wave functions of 



ref. [[331 have been used. The correlation function is taken from ref. |34| . 

The distorted wave functions in the final state are eigenf unctions of the optical potential 
of ref. f32[, which includes also the spin-orbit component. 

As explained in sect. 2, the nuclear current is the sum of a one-body part, whose contribu- 
tion is entirely due to the correlation function, and a two-body part, including contributions 
of the lowest-order diagrams with one-pion exchange, i.e. seagull, pion-in-flight and diagrams 
with intermediate isobar configurations. 

Calculations have been performed in the coplanar symmetrical kinematics, where the 
two nucleons are emitted in the scattering plane at equal energies and equal but opposite 
angles with respect to the beam direction, i.e. with ji = 72 and 712 = 71 + 72- In this 
kinematics, if E 1 is fixed and 71 is varied, it is possible to explore, for different scattering 
angles, all possible values of the recoil (j>b) or missing momentum (j>2m) distribution, where 
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P2m =PB = Q-Pl -P2 (23) 

If relative and CM motion are factorized and final-state interactions are neglected JE], p B is 
opposite to the total momentum P of the nucleon pair in the target. 



The choice of this kinematics is made at variance with that of ref. [[y| , where particular 
attention is drawn to the dynamics of dinucleons in the medium, which is studied through 
a kinematics where p B equals zero. The symmetrical kinematics is best suited to give 
information on the motion of the pair in different CM and relative angular momentum 
states and therefore to evaluate the effects of the one-body and two-body components of the 
nuclear current on these states, which is a goal of this paper. 

The calculated cross sections allow us to investigate a missing-momentum range up to 
about 250-300 MeV/c on both sides of the distributions. A wide range of photon energies, 
between 100 and 400 MeV, has been explored. 

The angular distributions corresponding to the five considered final states are given 
at E 7 =100 MeV in Figs. 1 and 2 and at E 7 =300 MeV in Figs. 3 and 4. In the figures 
separate contributions of the different terms of the nuclear current are compared with the 
cross sections produced by the sum of all of them. Moreover also the sum of the seagull 
and pion-in-flight currents and their sum with the one-body current are drawn, in order to 
better evaluate interference effects. The dependence on the photon energy in the considered 
range is given in Figs. 5 and 6, for the different contributions, at a fixed value of 71 near to 
the maximum of the distribution of each state. 

The shape of the angular distributions in Figs. 1-4 is determined by the combined effect 
of the CM wave function of the pair and of the different terms of the nuclear current. The 
main effect is given by the components of the CM wave function contributing to each state: 
+ and I2 have essentially an s-wave shape, if is due to a combination of p, d and s 
waves, while 2 + and 3 + are a combination of d and s waves at J5 7 = 100 MeV and have a 
d-wave shape at E y = 300 MeV, where the contribution of the A current is dominant. This 
analysis is consistent with the decomposition presented in Table I, for h.o. shell model wave 
functions, combined with the different values of the nuclear current matrix elements of the 
different relative states. 

The results in Figs. 1-4 indicate that the contribution of the one-body current and thus 
of SRC on the cross section is generally extremely small. It increases with the photon 
energy but is always overwhelmed by the two-body current. It is very sensitive to the choice 



of the correlation function: the strong correlation function of ref. [35] calculated with the 



hard-core OMY NN-interaction [3"6|, for instance, can enhance the contribution of the one- 



body current by about an order of magnitude. Even in this case, however, the calculated 
cross sections are dominated by two-body currents. The use in the present approach for pn 
knockout of a central and state-independent correlation function represents a very simple 
description of SRC. Moreover these correlation functions are not specifically calculated for a 
pn pair. Nevertheless, this simple treatment should be able to give a reasonable prediction 
of the role of SRC in the considered situations. Our numerical results indicate that (7,pn) 
cross sections are dominated by two-body currents and a more sophisticated treatment of 
correlations should not change this conclusion. This result is in agreement with that of 
refs. |T3|J2§. 



The seagull current is dominant, for all the considered transitions, at E 7 = 100 MeV and 
in general at low values of the photon energy, up to about 150 MeV. Its pure contribution 
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slightly decreases when the photon energy increases. The contribution of the pion-in-flight 
current is always much smaller than that of the seagull current and is only weakly dependent 
on the photon energy. It was found in refs. that a destructive interference is given by 

the sum of these two terms. This result is confirmed by the present calculations. However 
here this effect is generally small, but for the + state, i.e. the only state with T — 1, where 
a strong destructive interference is obtained. 

The isobar current shows, as expected, a strong dependence on the photon energy. Its 
contribution is negligible at low values of _E 7 , as in Figs. 1 and 2, but gains importance as 
the energy increases and becomes dominant above 200 MeV. For the four states with T = 
only the A-excitation current with its resonant propagator contributes. The resonance 
peak is clearly shown in the energy distributions of Figs. 5-6. On the contrary, only the 
A-deexcitation current contributes for the + (T = 1) state. Indeed for this state the 
A current does not show a resonant behaviour and above the pion-production threshold its 
contribution, although still dominant in the photon-energy range between 250 and 400 MeV, 
is smaller than for the other states. 

In Figs. 7-10 the photon asymmetries are shown in the same situations and in the same 
conditions as in Figs. 1-4 for the cross section. Different shapes are obtained for different 
states. The asymmetry is sensitive to the nuclear current, whose different effects on the two 
structure functions Wt and Wtt m Eq. are emphasized in the ratio. The differences 
between the results in Figs. 7 and 8, at E 1 = 100 MeV, and in Figs. 9 and 10, at E 1 = 300 
MeV, indicate the role of the different terms of the nuclear current at the two values of the 
photon energy. 

The asymmetry given by the contribution of the various terms of the current does not 
show a strong dependence on the photon energy. Like the cross section, it is not significantly 
affected by the one-body current. On the contrary, as it was already observed in the different 



model of ref. [15], it is very sensitive to all the components of the two-body current and to 



their interference. The pure contributions of the seagull and pion-in-flight currents generally 
give a positive asymmetry. In particular, for the seagull current it is positive and large (close 
to 1) for the + state and positive and small (around 0) for the 2 + state. A large effect is 
given by the interference between seagull and pion-in-flight terms, which generally produces 
a negative asymmetry. Only for the + state the asymmetry given by the sum of the two 
terms remains positive. This strong interference effect is essential to produce the final result 
at low values of E 1 , e.g. in Figs. 7 and 8 at E 1 = 100 Mev. The contribution of the pion-in- 
flight terms and of its interference with the seagull current is therefore much more important 
for the asymmetry than for the cross section. The asymmetry produced by the A current 
is generally negative for all the considered states. In the calculations, however, it presents 
anomalies and sharp peaks corresponding to minima of the cross section. At large values of 
the photon energy, e.g. in Figs. 9 and 10 at 300 MeV, the asymmetry, like the cross section, 
is dominated by the A current. However, for the + state, where only the deexcitation part 
of the isobar current contributes, the effect of the other terms of the two-body current is 
relatively a bit more important in determining the final result. 
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4. THE REACTION ie O(7,pp) 14 C 

For the reaction 16 0(7, pp) 14 C the theoretical approach of sect. 2 has been applied to the 
transitions to the three lowest discrete final states of the residual nucleus, i.e. the + ground 
state of 14 C, the 2+ state at 7.01 MeV and the 1+ state at 11.31 MeV. These are the states 
that have already been experimentally separated in the 16 0(e,e'pp) 14 C reaction |7]||[|. 

The two-nucleon overlap functions for each state are taken in the calculations from the 
spectral function of refs. fll3],|I2|] . Their explicit expressions, Eq. (Q), allow the separation 
of the different components of relative and CM motion. In Table II the main components 
and their amplitudes are given and compared with the corresponding amplitudes obtained, 
with h.o. wave functions, in the SM, where the states are described as two holes in the p 
shell: + is a (pi/2)~ 2 hole, while 2 + and 1 + are (P1/2P3/2) 1 holes. While the SM amplitudes 
are normalized to one, the spectral density amplitudes are normalized to the spectroscopic 
factor, i.e. to about 0.6. The table shows that in the spectral function the amplitudes 
for the 3 Pj relative states are more reduced than for the 1 Sq states with respect to the 
corresponding SM amplitudes. Other components, not shown in the table, contribute to the 
spectral density function and are included in the calculation, but their contribution is less 
than 1% HIH. 

The same conditions and the same ingredients used for the (7,pn) reaction in sect. 3 
have been adopted. For a pp pair, however, the two-body part of the nuclear current does 
not include the charge-exchange terms and contains only the non-exchange component of 
the A isobar current. 

In order to explore the missing-momentum distribution of each state, the same coplanar 
symmetrical kinematics as in sect. 3 has been considered, with the same range of photon 
energies and scattering angles. 

The angular distributions corresponding to the considered final states of 14 C are given in 
Figs. 11 and 12, at £ 7 = 100 MeV, and in Fig. 13, at E y = 300 MeV. In Fig. 11 the separate 
contributions of the one-body and A currents are compared with the cross sections given 
by the sum of the two terms. In Figs. 12 and 13 the contributions of the different relative 
states are drawn. 

The shape of the angular distributions is determined by the CM orbital angular mo- 
mentum of the knocked-out pair. The + state has essentially an s-wave shape; the p-wave 
component becomes more important at large values of the missing momentum, where the 
s-wave contribution becomes small. The 1 + state has a p-wave shape and the 2 + state is a 
combination of p and d waves. 

The A current is dominant in the cross section, for all the considered transitions and 
photon-energy values. Its contribution strongly increases with the energy, and above the 
pion-production threshold shows a pronounced resonant behaviour, with a maximum at the 
peak of the resonance. It is reduced in the 1 S state, where, as explained in sect. 2, the 
magnetic dipole NN «-» NA transition is forbidden. 

The contribution of the one-body current and therefore the effect of SRC is generally 
very small. SRC are accounted for in the calculations by defect functions which depend on 
the relative states. The results of Figs. 11-13 have been obtained with the defect functions 
calculated with the Bonn-A potential ||13|| . A different choice sensibly changes the contribu- 
tion of the one-boby current, but does only sligthly influence the final result, which is driven 



12 



by the A current. An example is shown in Fig. 14, where the cross sections obtained with 
the defect functions from Bonn-A and Reid Soft Core potentials [13| are compared, for the 
+ and 2 + states, at E 7 = 100 MeV where the contribution of the A current is smaller. At 
higher energies the one-body contribution increases and the results with Reid defect func- 
tions become larger than with Bonn-A. The increase with the energy of the contribution of 
the A current is however much larger and the final cross sections are even more dominated 
by the two-body current. Our results indicate that also (7,pp) cross sections are only slightly 
sensitive to SRC in the considered situations. 

In Fig. 15 the photon asymmetries are shown for the 1 + and 2 + states at E 1 = 100 
and 300 MeV. The contributions of the different relative states are shown with the same 
convention as in Figs. 12 and 13. For the + state, not drawn in the figure, the asymmetry 
is always negative and almost everywhere equal to —1. This result was already observed in 
ref. JE5 . 



The asymmetry, like the cross section, is dominated by the A current and does not show 
a strong dependence on the photon energy. Different results are obtained for different states. 
For 2 + the asymmetry is negative, while for 1 + it is large and positive, in agreement with 
the results of ref. [|15|], due to the large positive value (about 1) produced by the A current 
in the 3 Po state, which gives the dominant contribution to the final result. 

The asymmetry is not significantly affected by the one-body current. In Fig. 16 the effect 
of the different choices of the defect functions from Bonn-A and Reid potentials is shown for 
the 2 + state at E 1 = 100 MeV. It appears that the choice of the defect functions can have 
a significant but not substantial effect on the asymmetry, confirming that (7,NN) reactions 
are not particularly well suited for studying SRC. 



5. SUMMARY AND CONCLUSIONS 



In this paper exclusive 16 0(7,pn) 14 N and 16 0(7,pp) 14 C knockout reactions have been 
investigated, for transitions to the low-lying discrete final states of the residual nucleus. 

Exclusive (7,NN) reactions may represent a stringent test of reaction mechanisms and 
a promising tool for investigating the dynamics of bound pairs of nucleons in the nuclear 
medium. The presence of several discrete final states with well-defined angular momentum 
makes the 16 O nucleus an excellent target for this analysis. Measurements with sufficient 
energy resolution to separate individual final states and good statistical accuracy to deter- 
mine for each state the dependence of the cross sections on dynamical variables are needed. 
These data are now becoming available from new high-resolution experiments |23| . 

Numerical predictions of cross sections and photon asymmetries have been presented in 
the photon-energy range between 100 and 400 MeV. These results have been obtained in a 
direct knockout framework where both one-body and two-body currents are included. The 
contribution of the one-body current is entirely due to SRC; the two-body terms include 
contributions of the lowest-order Feynman diagrams with one-pion exchange, i.e. seagull, 
pion-in-flight and diagrams with intermediate isobar configurations, where a A is excited 
or a preformed A is deexcited after the absorption of the photon. Final-state interactions 
are taken into account by means of a phenomenological spin-dependent optical potential 
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describing the interaction of each of the outgoing nucleons with the residual nucleus. The 
two-nucleon overlap functions have been obtained for the 16 0(7,pp) 14 C reaction from a 
recent calculation of the spectral function JC!, where both long-range and short-range cor- 
relations are consistently treated. This approach was already applied to the 16 0(e, e'pp) 14 C 
reaction ||12|| . The two-nucleon overlap functions for the 16 0(7, pn) N reaction are given by 
the product of a coupled and fully antisymmetrized SM pair function and a Jastrow type 
correlation function, which incorporates SRC, since a calculation of the spectral function for 
a pn pair is not available. Although presumably too simplistic, this prescription should be 
able to give reasonable predictions of the role of SRC in (7,pn) reactions. 

The angular distributions of the cross sections for transitions to different states have 
different shapes, essentially determined by the components of the CM orbital angular mo- 
mentum of the initial pair. Also the asymmetry is sensitive to the behaviour of the different 
components of the pair wave functions. These results confirm that exclusive reactions are 
particularly well suited to check the relevance of the direct knockout mechanism and to 
explore the conditions of different pairs of nucleons in the nuclear medium. 

Cross sections and asymmetries, for both (7,pn) and (7,pp) reactions, are only slightly 
affected by SRC and are dominated, in all the considered situations, by two-body currents. 
Therefore photon-induced two-nucleon knockout reactions do not seem particularly well 
suited to investigate SRC, but might give interesting information on the various terms of 
the nuclear current and on their behaviour in different kinematic conditions and in different 
states. 

Better information on correlations for a pn pair might be obtained from the (e,e'pn) 
reactions, where also the longitudinal component of the nuclear response contributes. In this 
case, as it was already observed in (e,e'pp) [l^] , a suitable choice of kinematics in exclusive 
experiments might allow the separation of effects due to SRC and two-body currents. For this 
analysis, however, a more refined treatment of correlations would be desirable, such as that 
applied here to pp pairs through the spectral function. Moreover, also tensor correlations 
might be investigated for nucleon pairs in a T = state. 

Our results indicate that in the 16 0(7,pn) 14 N reaction the seagull current is dominant 
at lower values of the photon energy, up to about 150 MeV, then the A current gains 
importance. The contribution of the pion-in-flight current is always much smaller than 
that of the seagull current and is only weakly dependent on the photon energy. In the cross 
sections it generally produces a destructive interference with the seagull currents. This effect 
is however small, but for the + (T = 1) state. Interference effects are better emphasized 
in the photon asymmetry, where the interference between seagull and pion-in flight terms 
turns the asymmetry from positive, for both separate contributions, to negative, but for the 
+ state. This effect is essential at low values of the photon energy, while at higher energies 
both cross sections and asymmetries are dominated by the A current. The A current is 
dominant also in the 16 0(7,pp) 14 C reaction, for all the considered transitions and over all 
the considered range of photon energy. Its contribution strongly increases with the energy, 
but plays the main role, both on cross sections and asymmetries, even at 100 MeV. 

The results of our investigation indicate that a good description of the isobar current 
is indispensable in the theoretical treatment. On the other hand, (7,NN) reactions might 
give deeper insight into the behaviour of the A in the nuclear medium. Moreover, also the 
influence of two-body currents due to heavier meson exchange could deserve further inves- 



14 



tigation. 

We are grateful to K. Allaart and P. Grabmayr for useful comments and fruitful discus- 
sions. 
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TABLE I. Proton-neutron removal amplitudes from 16 for states of 14 N in a h.o. approxi- 
mation and for different relative 2S+ 'lj states. L is the CM angular momentum. 
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TABLE II. Two-proton removal amplitudes from 16 for states of 14 C in a h.o. approximation 
and for different relative 2S+ 'lj states. L is the CM angular momentum. The amplitudes are 
compared with the corresponding coefficients c 1 in Eq. (||), obtained from the calculation of the 
spectral function of ref. |L| . 
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FIGURES 



FIG. 1. The differential cross section of the reaction 16 0(7, pn) 14 N, as a function of the angle 
71, in coplanar symmetrical kinematics at = 100 MeV, for the transitions to the low-lying states 
in 14 N: if (E 2m = 22.96 MeV), 0+ (E 2m = 25.27 MeV), if (E 2m = 26.91 MeV), 2+ (E 2m = 29.99 
MeV), 3 + (E 2m = 34.01 MeV). The optical potential is taken from ref. |3^], the single-particle 
wave functions from ref. [^] and the correlation function from ref. |34|| . Separate contributions of 
the one-body current (1-B) and of its sum with the two-body seagull (SEA) and pion-in-flight (ir) 
currents are drawn. The solid lines give the total cross section, where also the A current (A) is 
added. 

FIG. 2. The differential cross section of the reaction 16 0(7, pn) 14 N, as a function of the angle 71, 
in the same kinematics and for the same transitions as in Fig. 1. Optical potential, single-particle 
wave functions and correlation function as in Fig. 1. The solid lines are the same as in Fig. 1. 
Separate contributions of the A, seagull, pion-in-flight currents and of the sum of seagull and 
pion-in-flight currents are drawn. 

FIG. 3. The differential cross section of the reaction 16 0(7, pn) 14 N, as a function of the angle 
71, in coplanar symmetrical kinematics at E~ = 300 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 1. 

FIG. 4. The differential cross section of the reaction 16 0(7, pn) 14 N, as a function of the angle 
71, in coplanar symmetrical kinematics at E~ = 300 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 2. 

FIG. 5. The differential cross section of the reaction 16 0(7, pn) 14 N, as a function of the photon 
energy, for a fixed value of the angle 71, in coplanar symmetrical kinematics, for the same transitions 
as in Fig. 1. Optical potential, single-particle wave functions and correlation function as in Fig. 1. 
Line convention as in Fig. 1. 

FIG. 6. The differential cross section of the reaction 16 0(7, pn) 14 N, as a function of the pho- 
ton energy, for a fixed value of the angle 71, in coplanar symmetrical kinematics, for the same 
transitions, under the same conditions and with same line convention as in Fig. 2. 

FIG. 7. The photon asymmetry of the reaction 16 0(7, pn) 14 N, as a function of the angle 71, in 
the coplanar symmetrical kinematics at EL = 100 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 1. 

FIG. 8. The photon asymmetry of the reaction 16 0(7, pn) 14 N, as a function of the angle 71, in 
the coplanar symmetrical kinematics at E~ = 100 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 2. 
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FIG. 9. The photon asymmetry of the reaction 16 0(7, pn) 14 N, as a function of the angle 71, in 
the coplanar symmetrical kinematics at EL = 300 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 3. 

FIG. 10. The photon asymmetry of the reaction 16 0(7, pn) 14 N, as a function of the angle 71, in 
the coplanar symmetrical kinematics at EL_ = 300 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 4. 

FIG. 11. The differential cross section of the reaction 16 0(7, pp) 14 C, as a function of the angle 
71, in coplanar symmetrical kinematics at EL = 100 MeV, for the transitions to the low-lying 
states in 14 C: 0^ (E 2m = 22.33 MeV), 2+ (E 2m = 30.00 MeV), 1+ (E 2m = 33.64 MeV). The 
defect functions for the Bonn-A potential and the optical potential of ref. |32| are used. Separate 
contributions of the one-body and two-body A currents are drawn. The solid lines give the total 
cross sections, produced by the sum of the one-body and A currents. 

FIG. 12. The differential cross section of the reaction 16 0(7, pp) 14 C, as a function of the angle 
71, in the same kinematics and for the same transitions as in Fig. 11. Defect functions and optical 
potential as in Fig. 11. The solid lines are the same as in Fig. 11. Separate contributions of different 
partial waves of relative motion are drawn. 

FIG. 13. The differential cross section of the reaction 16 0(7, pp) 14 C, as a function of the angle 
71, in coplanar symmetrical kinematics at E 7 = 300 MeV, for the same transitions, under the same 
conditions and with same line convention as in Fig. 12. 

FIG. 14. The differential cross section of the reaction 16 0(7, pp) 14 C, as a function of the angle 
71, in coplanar symmetrical kinematics at E 7 = 100 MeV, for the transitions to the + and 2 + 
states in 14 C. The solid and dashed lines are calculated with the defect functions of the Bonn-A 
and Reid potentials, respectively; the dotted and dot-dashed lines give the corresponding separate 
contributions of the one-body current. The optical potential is taken from ref. |3^] 

FIG. 15. The photon asymmetry of the reaction 16 0(7, pp) 14 C, as a function of the angle 71, for 
the transitions the 2 + and 1 + states in 14 C in the same kinematics as in Fig. 12. Defect functions, 
optical potential and line convention as in Fig. 12. 

FIG. 16. The photon asymmetry of the reaction 16 0(7, pp) 14 C, as a function of the angle 71, 
in coplanar symmetrical kinematics at E 7 = 100 MeV, for the transition to the 2 + state in 1 C. 
Optical potential and line convention as in Fig. 14. 
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